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ABSTRACT 

We present a new technique for detecting scattered starlight from transiting, cfose-orbiting ex- 
trasolar giant planets (CEGPs) that has the virtues of simplicity, robustness, linearity, and model- 
independence. Given a series of stellar spectra obtained over various phases of the planetary orbit, 
the goal is to measure the strength of the component scattered by the planet relative to the com- 
ponent coming directly from the star. We use two complementary strategies, both of which rely on 
the predictable Doppler shifts of both components and on combining the results from many spectral 
lines and many exposures. In the first strategy, we identify segments of the stellar spectrum that are 
free of direct absorption lines and add them after Doppler-shifting into the planetary frame. In the 
second strategy, we compare the distribution of equivalent-width ratios of the scattered and direct 
components. Both strategies are calibrated with a "null test" in which scrambled Doppler shifts are 
applied to the spectral segments. As an illustrative test case, we apply our technique to spectra of HD 
209458 taken when the planet was near opposition (with orbital phases ranging from 11 to 34°, where 
0° is at opposition), finding that the planet-to-star flux ratio is (1.4 ± 2.9) x 10""' in the wavelength 
range 554—681 nm. This corresponds to a geometric albedo of 0.8 ± 1.6, assuming the phase function 
of a Lambert sphere. Although the result is not statistically significant, the achieved sensitivity and 
relatively small volume of data upon which it is based are very encouraging for future ground-based 
spectroscopic studies of scattered light from transiting CEGP systems. 
Subject headings: planetary systems — stars: individual — (HD 209458) — techniques: spectroscopic 



1. INTRODUCTION 

The discovery of the first exoplanet orbiting an or- 
dinar y main sequence star, 51 Peg, a dozen years 
ago (jMavor fc Oueloj Il995( l also revealed the first 
example of a largely unexpected but common class 
of astronomical objects, namely close-orbiting extra- 
solar giant planets (CEGPs). The discovery of 

the first transi ting exoplan et sys tem, HD 209458b 

also a CEGP (iHenrv et al.l l2000t iMazeh et al.] 120001 : 
ICharbonneau et al.ll2000l : iBrown et al.ll2001[ ). a few years 
later allowed study of these enigmatic objects to move be- 
yond orbital parameters and minimum masses. Although 
radial velocity, transit and microlensing techniques have 
made the study of exoplanets one of the most active 



and exciting areas of astrophysics in recent years and 
have produced a wealth of additional major discover- 
ies and import ant theoretical problems ( Perrvman '2000|; 
IChauvinll2007| ). the effort to understand CEGPs remains 
central to the field. 

The light of the primary star scattered from CEGPs 
reflects the nature (composition, clouds, energy balance, 
structure, motions et cetera ) of their atmospheres (e.g., 
Seager fc Sasselov|[T998: Marlev et al.iil999: Seager et all 



20001: ISudarskv et al.ll2000L 120031: iBurrows et al.ll2005li : 
its direct detection would thus test and discriminate 
among emerging models of CEGPs. 

There have been previous attempts to detect scat- 
tered starlight from CEGPs against the much brighter 
direct light from the primary via high-resolution spec- 
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troscopic searches for the scattered, and hence Doppler- 
shifted, stellar absorp tion lines (jCharbonneau et all 
Il999t iLeigh et al.1 l2003l and references therein). The 
problem is very difficult because the expected contribu- 
tion of the scattered light to the total light is of order 
10~^ or less for known CEGPs. Due to both the very 
small amplitude of the expected signal and the unknown 
orbital inclination angle i for non-transiting systems, the 
methods so far employed rely on model-dependent as- 
sumptions to estimate prior distributions for planetary 
radii, orbital inclination s and the form of the phase func- 
tion (jLeigh et al.ll2003f l. The advantage of such an ap- 
proach is that it can be applied to the brightest (but non- 
transiting) CEGP systems that thus have the least pho- 
ton noise. In other words and as is usually the case, rely- 
ing on a prior model allows one to maximize purely sta- 
tistical signal-to-noise (S/N) but at the price of increased 
systemat ic uncertainty. Two previously published meth- 
ods (e.g. lCharbonneau et al.|[l999l : iCameron et aLlll999D 
depend on a x^ minimization applied to the whole stel- 
lar spectrum, which is not an optimally sensitive method 
for this problem because it dilutes the detectability of the 
scattered lines (which carry the signal of interest) by con- 
sidering the quality of the fit to the continuum regions. 

The simplest example of the advantage of studying 
transiting CEGP systems is that measurement of the 
scattered light fraction allows direct determination of 
the albedo in isolation, and not merely in combination 
with the planet's radius squared (which must be esti- 
mated from a theoretical model for non-transiting sys- 
tems). The major advantage of studying non-transiting 
systems is that some of them are nearly two orders of 
magnitude brighter than the brightest transiting system, 
thus yielding far smaller photon statistics noise. 

Here we propose a largely model-independent tech- 
nique to detect the scattered stellar absorption lines 
in transiting CEGP systems. In this method, no as- 
sumptions concerning the form of the phase function, 
the planet's radius, the orbital inclination et cetera are 
needed to measure the planet-to-star flux ratio and 
albedo. Rather, only physical parameters of the system 
measured via transit observations are required. More- 
over, the technique maximizes sensitivity by concentrat- 
ing only on those portions of the spectrum in which scat- 
tered stellar lines are known to be present (at some am- 
plitude to be measured). As a practical "worked exam- 
ple" we apply this technique to Subaru HDS spectra of 
HD 209458. The required input physical parameters of 
the system and their relation to the expected scattered 
light signal are briefly reviewed in §2. In §3, we describe 
the technique in detail, and apply it to HD 209458b, us- 
ing high precision optical spectra acquired on UT 2002 
October 27, with the Subaru Hig h Dispersion Spectro- 
graph (HDS; iNoguchi et al.l[200^ . Results are summa- 
rized in §4. Prospects for future studies are then exam- 
ined based on these results. 

2. INPUT PARAMETERS 

Our technique utilizes known physical properties of the 
transiting CEGP system based on prior radial velocity 
and transit studies. 

The planet's velocity re lative to the star at orbita l 
phase (j) is given by (see, e.g. lColfier Cameron et al.ll2002l ) 

Vp{(j)) ^Kp-sm2TT(j), (1) 



where 4> = {t — To)/Porb is the orbital phase at time t, Tq 
is the transit epoch, and Port is the orbital period. Its 
apparent radial-velocity amplitude Kp about the mass 
center of the system is 



Kp^ 



27ra sini 



^orb -L 



q 



(2) 



where a is the orbital distance, i is the inclination and 
q = Mp/Mi, is the mass ratio. The values of a and i can 
be measured from radial- velocity and transit-photometry 
data, respectively. 

If the planet is tidally locked, as expected for nearly 
all CEGPs, there will be no rotational broadening. Then 
the scattered spectra will be the Doppler-shifted stellar 
spectra, scaled by a factor of the pla net-to-star flux ratio 
e, which is usually decomp osed as (jCharbonneau et all 
Il999t ICameron et al.|[l99^ 



i 



e(a,A) =_p(A)5(a,A) 



i?p 



(3) 



where p{X) is the geometric albedo, and g{a, A) is the 
phase function, with a being the phase angle which varies 
as cosa = — sini cos27r0. The value of Rp is also mea- 
sured from transit-photometry data. 

3. METHOD AND APPLICATION TO HD 209458B 
3.1. Overview 

The basic idea behind the method is that for any ab- 
sorption line in the primary star's spectrum we can pre- 
cisely predict the wavelength of each line in the light 
scattered from the planet in each exposure; it is simply 
offset from the position of the stellar line by the Doppler 
shift due to the planet's orbital motion relative to the pri- 
mary, as projected onto the line of sight to the observer 
at the time of that exposure. The strength of the scat- 
tered line relative to the line in the primary's spectrum is 
just the fractional contribution of the planet's scattered 
light to the total light at that wavelength; in principle, 
this line strength gives a measure of the planet's albedo 
at that wavelength and phase angle. 

Of course, in practice, for known CEGP systems this 
fraction is expected to be of order 10~^ to 10~^ and any 
individual line is thus completely undetectable even in 
very high S/N spectra. Our technique is then simply to 
superimpose and combine pieces of the spectra at the 
expected wavelengths of the scattered lines over many 
lines and exposures to achieve sufficient S/N to allow a 
statistically significant detection. This can be done either 
1) by simply summing small fragments of the spectra 
centered on these predicted locations or 2) by measuring 
the ratio of scattered absorption line equivalent-widths 
to the corresponding equivalent widths in the primary 
and then combining the ratios over lines and exposures. 
We employ both variations. 

In an attempt to detect tiny fractional signals against 
a noisy background, it is obviously crucial to know the 
background rather precisely and to characterize its noise 
properties accurately. We do both by an "internal" pro- 
cedure using the same spectra and same wavelength re- 
gions in which we search for the signal. Specifically, 
in addition to combining the spectral fragments defined 
by the expected wavelength positions of the scattered 
lines, we also carry out the same analysis but with the 
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Doppler shift assigned to each exposure "shuffled" ran- 
domly among the exposures. In this way we can produce 
a large number of reference spectra made up of the same 
spectral fragments taken from the same set of exposures 
but each lacking any scattered line contribution. By com- 
bining this ensemble of reference spectra we obtain an ex- 
cellent estimator of both the background against which 
we must distinguish the scattered light signal and of all 
of its noise properties (including, for example, the collec- 
tive effects of weak stellar lines that cannot be detected 
individually) . 

Finally, in order to calibrate and test the method, it 
is very useful to have mock or simulated data containing 
an artificially injected, and thus precisely known, signal 
and having realistic noise properties. This we also man- 
ufacture internally. The individual spectra of the sys- 
tem are shifted to the planet's rest frame at each expo- 
sure, and then multiplied by some assumed planet-to-star 
fl-ux ratio e to produce a set of artificial planet scattered 
light spectra. These artificial spectra are then added 
back into the individual exposures. For example, if the 
original spectra contain some scattered light signal with 
a planet-to-star ratio eo, then these artificially injected 
spectra should have scattered light with a planet-to-star 
ratio ei = (eo + e)- We can analyze the artificially in- 
jected spectra to obtain ei and then subtract eg from it 
to recover e. Successful, or failed, recovery of this known 
artificial signal by application of the same technique pro- 
vides a strong test of the reliability and sensitivity of the 
method. 

3.2. Data Preprocessing 

As an illustrative case study, we analyze 33 Subaru 
HDS spectra of the system HD 209458, taken on UT 
2002 October 27, when the planet was just out of the sec- 
ondary eclipse, with the orbital phase 11.0° < a < 33.9°. 
A detailed description of the observations and data re- 
d uction to obtain on e-di mensional spectra c an be found 
in lWinn et all (|2004D and lNarita et al.l(|2005f) . The wave- 
length coverage of the spectra analyzed is from 554 nm 
to 681 nm, with a spectral resolution of i? « 45000. The 
typical exposure time was 500 sec, and the S/N per pixel 
was ~ 350. 

The reasons for adopting exposure phases in the vicin- 
ity of the secondary eclipse, but not extremely close to it, 
are two fold. First, the strong back scattering of aerosol 
clouds or grains, if any, could dramatically brighten the 
planet. Second, we should not observe too close to ei- 
ther the primary transit or the secondary transit, be- 
cause during those events the planet's motion is nearly 
in the plane of the sky and its relative radial velocity vari- 
ation is small, causing the scattered absorption lines to 
be blended with the direct absorption lines. For the HD 
209458 data analyzed here, the velocity offsets between 
the stellar lines and their scattered counterparts is in the 
range of 27 — 80 km s~^, sufficient to avoid blending in 
all cases. 

The one-dimensional spectra are first normalized. We 
fit an ll*'*-order-spline3 function to obtain the blaze func- 
tion of each order separately, and then merge the blaze 
functions of all the orders. The un-normalized spectra of 
each order arc also merged in the same way and divided 
by the merged blaze function. Spectra are re-sampled 
and dispersion-corrected to logarithmically spaced bins 



TABLE 1 

Physical para meters of the system HD209458 from 

IKnutson et al.I (J2007I '). 



Parameter 



Value and references 



Mass of the star: Ah (Mq) 
Radius of the star: _R* (Rq) 
Mass of the planet: Mp (Mj) 
Radius of the planet: Rp (Rj) 
Orbital inclination: i (degrees) 
Orbital period: Porb (days) 
Transit epoch: Tq (HJD) 



, ,(,,+0.066 

1.1U1_Q Qg2 

1 19C-+0.020 

l.lZO_Q Q23 

0.64±0.06 

1 Q9n+0024 
l.ozu_Q Q25 

86 929+00°'' 

3.52474859±0.00000038 

2,452,826.628521±0.000087 



after merger. All exposures are normalized separately, as 
their blaze functions vary due to inst rumental variations 
(jWinn et al.ll2004l : iNarita et al.ll2005f ). We note that the 
method to be discussed in Section !??^ is also applicable to 
the un-normalized spectra. The normalization however 
simplifies the EW test to be described in Section 13.41 
The reason we use an ll*^-order-spline3 function for the 
fitting is to remove the large-scale profiles of the spectra 
without diluting the weak lines containing the scattered 
light signal. 

We adopt the physical parame ters of the system HD 
209458 from lKnutson et al.l(|2007D . as summarized in Ta- 
ble m along with the Modified Julian Day (M JD) at the 
central time of each exposure, to predict Vp ((/>). 

3.3. Strategy I: Summing Scattered Light Lines 

In this section, we search for the scattered light by 
aligning and combining the imprinted signals over a large 
set of stellar absorption lines as well as over exposures. 
The scattered signals are aligned by Doppler-shift cor- 
recting each exposure according to the corresponding 
Vp{4)). We also make a set of reference spectra by com- 
bining the same set of exposures, which, on the other 
hand, are Doppler-shift corrected with the Vp{(j)) of an- 
other exposure. Therefore, extra absorptions around the 
predicted scattered-light centers would indicate a detec- 
tion. 

3.3.1. Details of the Strategy- 1 Procedure 

First, 33 pre-processed spectra Sy ^^ are Doppler-shift 
corrected according to their predicted corresponding ra- 
dial velocities. We use a logarithmic wavelength scale, 
and the super script "0" to refer to the star's rest frame 
and "Z" for the planet's. Then the shifted spectra, S^ )^, 
are corrected to have the same dispersion as that of S^ ^^ . 
Here and throughout this study, spectra are always cor- 
rected to this standard dispersion after being shifted. 

Second, a set of suitable stellar absorption lines were 
selected according to the following criteria: 1 - The line's 
full-width at half-maximum (FWHM) must be larger 
than the threshold value FWHMt/j = 0.15A, a value cho- 
sen by experimentation to maximize sensitivity. (The 
final result is, however, rather insensitive to the exact 
value of FWHMt/i.) 2 - The line's blue^-ends extending 
to a wavelength range, which corresponds to the largest 

^ Due to the planet's phases, its scattered spectra in all the 33 
exposures we used in this study would be blue-shifted relative to 
the stellar spectra. 
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Fig. 1. — The 123 selected stellar absorption lines. Wavelength 
versus FWHM (top), and the distribution of FWHM {bottom) . 

radial velocity among all the exposures, needs to be rel- 
atively free of strong stellar lines. This is required to 
avoid contamination of any detection by strong stellar 
lines. The wavelengths of selected lines and the distribu- 
tion of their FWHMs are displayed in Figure [TJ 

Third, small fragments (AA ~ 1.7 A each^) of the 
spectra are averaged together, centering on the predicted 
scattered features and over all the selected lines and ex- 
posures. The combined normalized counts as a function 
of wavelength difference (from the predicted line centers) 
can be written as 



C^{k) 



1 



JVe N, 



NeNl 



^^5^=[log(Afe,,0-log(Ao,eO], (4) 



e=l 1 = 1 



where k = —50, —49, ...,49, 50, and represents the fc-th 
pixel on the bluer (redder) side of the scattered-light 
center (fc = 0), if fc is negative (positive); log(Afc+i_ez) — 
log(Afc,ei) corresponds to dX ~ 0.017 A, i.e., the wave- 
length pixel size on the detector at the center of the 
spectral coverage'^; N^ — 33 is the number of exposures, 
whereas Ni = 123 is the number of selected stellar ab- 
sorption lines. Ze stands for the "redshift" of the e*^ 
exposure, and S^' denotes the spectrum that has been 
Doppler-shift corrected according to the radial velocity of 
the e*^ exposure. Wavelength solutions of different expo- 
sures agree with each other to within ^0.03 A, therefore 
the central wavelength of each stellar line, Ao,ei, is a func- 
tion of both line and exposure. We account for the small 
variations in wavelength solutions by localizing each line 
in each exposure separately. In other words, C^{k) is 
simply the average scattered spectral line generated by 
summing, after continuum normalization, over all of the 
selected lines (123 for our HD 209458 data) and all of 
the exposures (33 for our HD 209458 data), a total of 

^ A A PS 1.7 A was chosen to be big enough to cover several 
typical stellar line widths. 

^ The actual spectral resolution at this point is fs 0.14 A, equiv- 
alent to just about 8 detector pixels. 



Ne X Ni separate spectral fragments (4059 in the case of 
the present data). 

Fourth, we take the mean (as well as the median) 
values calculated from the same 33 pre-processed spec- 
tra, Doppler-shift corrected using N sets of scrambled 
redshifts (hereafter, reference spectra, for short), as the 
standard to compare C^ {k) with. This spectrum can be 
regarded as the "continuum" level, i.e., relatively free of 
scattered signals, and is given by 



1 ^ 



(5) 



^N^ 



^N, 



Where C^.(fc) = jr7Nl^e^iU=iS^--^[\og{X,^,i) ~ 

log(Ao,ei)]- {Z'ne} '^^ ^ permutation of {Ze}, and there 
are A^ = 1, 000 sets of random permutations. We note it 
is essential that {Z,'j g} is a permutation of {Z^}, rather 
than purely random redshifts, because using a permuta- 
tion of {Ze} ensures the exact same fragments of stellar 
spectra go into C^"(fc), as those in C^(fc). 

Fifth and finally, we define the combined scattered sig- 
nal as the residual given by 



AC(fc) EE C^(fc) - (C^ (k)). 



(6) 



Given scattered light, AC{k) will be negative around 
fc = 0, with a profile similar to the combined stellar lines 
scaled by e(a. A). Strictly speaking, (C^ (fc)) is affected 
to some degree by the scattered-light signal, even though 
the redshifts have been scrambled. This is because oc- 
casionally the randomly assigned redshift will happen to 
be or close to the appropriate redshift of the exposure. 
Therefore (C^ (fc)) serves as a lower limit for the true 
continuum level, and AC(fc) yields a conservative value. 

3.3.2. Results for Strategy I 

Figure [2] shows the normalized count C(fc) and 
the residual AC(fc) obtained from analyzing the Sub- 
aru/HDS spectra of HD 209458b. The overall shape 
of C(fc) is a result of the combination of multiple stel- 
lar spectrum segments. The significance levels plotted 
in (A), as well as the 1-a errors shown in (B) are cal- 
culated from the variance among the reference spectra. 
AC(fc) along with the 1-ct errors in (C) and (D) are from 
the jackknife- and bootstrap- re-samplings (re-sampled 
over both lines and exposures) of the data along with 
the corresponding reference spectra, respectively. The 
re-sampling tests produced similar errors to those solely 
estimated from the variations in reference spectra. To 
test the algorithms, we also analyzed a mock/simulated 
spectrum, obtained by combining the same set of stellar 
lines scaled by a factor of 1.2 x lO"**. This simulated 
signal indicates the expected scattered light assuming a 
planet-to-star fiux ratio of e = 1.2 x 10~*. The com- 
parison of the observations and simulated signal shown 
in Figure [2] is intriguing in that the data are consistent 
with a low statistical significance detection of the scat- 
tered light, with e > 10"** at the l-tr level, as the data 
follow the simulated signal tantalizingly well, especially 
around the predicted line center where the signal should 
be the strongest. 

Figure [3] also displays the normalized residual counts 
AC(fc), but in this figure the actual spectral resolution. 
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Fig. 2. — The normalized count C{k) and the residual AC(fc) as a function of wavelength difference k, relative to the predicted scattered 
line's center, aligned and average-combined over a large set of stellar lines and over exposures. Each exposure has been Doppler-shift 
corrected to align the scattered signals, according to the planet's relative radial velocity to the star. (A): the red curve shows C^{k), 
whereas the black curves represent the median (solid; also same as the mean) and 68% (dotted), 95% (dashed), and 99% (dash-dotted) 

significance levels, obtained from 1,000 sets of scrambling-shifted spectra ({C "(fc)}); (B), (C) and (D); the black filled-circles show the 
residual AC(k), which is the difference between C^ (k) (red) and (C^ (k)) (black-solid). Results from the simulated scattered light by 
combining the stellar lines scaled by a factor of 1.2 X 10~* are shown as blue-dashed curves. The error bars in (B) correspond to 1 a- 
errors obtained from the significance levels shown in (A), whereas (C) and (D) present the results of 4,059 sets of jackknife- and bootstrap- 
re-samplings. In all the panels, one pixel corresponds to ~ 0.017A, the approximate size of a detector pixel, about one eighth of the actual 
spectral resolution. The systematic decrease in AC(fc) at more than 12 pixels to the red or blue of the line center is due to the influence 
of stellar absorption lines beyond the edges of the region chosen to bo clear of such interference. 
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Fig. 3. — The normalized residual counts AC'(k) as a function of wavelength difference k, as in panels (B) and (D) of the previous figure, 
but in this case the data points have been rebinned by a factor of 8. This corresponds to the true spectral resolution of the data, rather 
than to the detector pixel size. The rebinnig is arranged such that the predicted line center position at fc = is in the middle of the central 
point. The dashed line shows the expected shape of the composite scattered light line assuming a 1.2 X 10"* scattered light fraction. 
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about 8 detector pixels, is used. There are thus only three 
points across the expected composite scattered light ab- 
sorption line, but the error bars are correspondingly re- 
duced, and the points are not strongly correlated by the 
finite spectral resolution as they are in Figure [2l Since 
the two figures simply display the same data at two differ- 
ent wavelength resolutions, they are both naturally con- 
sistent with the same scattered-light fraction of 1.2 x 10"'* 
and the same line shape, shown by the blue-dashed line. 

3.4. Strategy II: Equivalent Width Ratio Distributions 

In this section, we further suggest a method to mea- 
sure the planet-to-star flux ratio e by examining the dis- 
tribution of the EW ratios between the scattered and the 
stellar absorption lines, measured for a large set of lines 
and exposures. The FWHMs of stellar lines displayed in 
Figure [1] are used to set the wavelength range over which 
the EWs of both stellar and scattered lines are evalu- 
ated. Since this method effectively concentrates all of 
the "missing flux" signal into a single number with appro- 
priate weighting of the stronger and weaker features, we 
would expect better sensitivity performance than Strat- 
egy I. The cost, of course, is loss of all information about 
the shape of the composite shape of the scattered lines 
and thus of a useful check that the measured feature has 
the expected character of scattered light. The two ap- 
proaches are thus complementary. 

3.4.1. Details of the Strategy- II Procedure 

We adopt a dimensionless EW, which is convenient for 
a logarithmically linear wavelength scale, given by 



For the scattered signals, W^ 



A ■ EJi 



w = — 



dX 



1- 



_Sx_ 
S{0) 



= lnlO- 



1- 



Sloe- 



S{0) 



rflogA, (7) 



where Aq is the central wavelength of a stellar absorption 
line and S{0) denotes its local continuum. In logarith- 
mically linear wavelength scale, W can be discretized as 



w = a-Y: 



1- 



3(0) 



(8) 



where A = InlO • AlogA is a constant, and the summa- 
tion goes over all pixels inside a wavelength range of twice 
that of the line's FWHM and centered at Aq. Since dif- 
ferent stellar lines have different FWHMs as shown in 
Figure [U a variable number of pixels are included in the 
EW integrations by making use of the line-strength in- 
formation, whereas in Strategy I, the same number of 
pixels are combined for all lines. 
The EW of the stellar lines, M^*; = A • ^Jl - 

Si Xi ci/^eii^)]^ i^ calculated for each hue (denoted by 
/) in each exposure (denoted by e), using the stellar 
FWHMs measured in the initial exposure, which are 
shown in Figure [TJ To check for FWHM variability, we 
measure the FWHMs in both the initial and final expo- 
sures and find that they agreed within 0.1 per cent. 

For each line, the local continuum 5*^ ;(0) is estimated 
by averaging the median of the counts in two spectral 
intervals on either side of the line, with a size chosen to 
be neither too large to reflect the local continuum level 
nor too small for accurate statistics. We use a range of 
50 pixels, corresponding to ^1 A. The results have been 
tested to be independent of the adopted size of the range 
over which the continuum levels are estimated. 



■^logAi ef/'^f; (*^)]' where the super script Z^ represents 
the spectrum S has been Doppler-shift corrected with the 
Doppler shift of the e*'' exposure. We use the FWHMs 
of the stellar lines, assuming that the scattered light is 
a replica of starlight, scaled by a factor of e. The local 
continuum S^ f (0) is taken as the median of the line's 
blue end, since the red end often overlaps the associ- 
ated strong stellar lines. The EW ratio is then given by 
R^l = W^l/Wli- The EW-ratio measurements of Ni 
lines in N^ exposures constitute a statistical sample of 
{i?g;}, which has Ni x N^ data points. 

z' 
We also calculate {R^ ;'} for TV sets of reference spectra, 

to construct a standard comparison pool. Both mean 

and median are adopted as statistical measures, where 

M{R^) denotes taking either the mean or the median of 

{i?f;}. Finally, we define the residual EW ratio, given 

by ' 



zs 



AAf (i?) = M{R^) - {MiR-^ )}, 

1 ^ 
where (Af (i?^ )) = - ^ A/(i?^"), (9) 

n=l 

is the comparison standard obtained from N sets of 
scrambling-shifted reference spectra. Given null (by con- 
struction) reference-spectrum signal, AM{R) will be pos- 
itive and give a direct estimate of e. 

3.4.2. Results for Strategy II 

Figure m shows AM{R) for both the median and the 
mean statistics, obtained from the HD 209458b data, 
simulated scattered signals, reference spectra, and boot- 
strap re-samplings. We have also generated jackknife 
re-samplings, which produced similar results to those of 
the bootstrap re-samplings, and both re-sampling analy- 
sis yield consistent errors with those estimated from the 
variation of reference spectra. To test the linearity of 
the algorithm, we injected mock/simulated scattered sig- 
nals, assuming e = 10"'* and 2 x 10"'*, respectively. The 
injected signals are recovered linearly through the data 
pipeline as expected, since our algorithm only involves 
linear transformations. These EW-ratio results suggest 
e = (1.4 ± 2.9) x lO""* ((3.6 ± 2.0) x 10"'*), according to 
the median (mean) statistics. The mean is more sensitive 
to outliers which can be caused by stellar-line overlaps 
or inaccurate estimation of the local continuum levels, or 
by real intrinsic variation of the albedo with wavelength 
or phase, whereas the median should be more robust and 
representative of a global effective albedo over the whole 
wavelength and phase ranges. 

4. SUMMARY AND DISCUSSION 

We have developed a simple, robust, linear, and model- 
independent technique to detect scattered starlight from 
transiting CEGPs, by measuring scattered spectral fea- 
tures at predictable Doppler shifts from the stellar lines 
due to planets' orbital motions. There are two main vari- 
ations of this technique. First, spectrum fragments have 
been aligned and combined over a large set of stellar ab- 
sorption lines as well as over many exposures taken near 
the secondary eclipse. Second, the planet-to-star flux ra- 
tio has been determined by examining the distribution 
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Fig. 4. — The EW ratio AM{R). (A) and (C) show the mean results whereas (B) and (D) present the median. In (A) and (B), the 
histogram displays the normalized distribution of results from 1,000 sets of scrambling-shifted spectra. The dotted line represents the 
averages of the distributions. As shown by the legend in (B), the other three vertical lines show results calculated from the data [thick red- 
solid), the simulated scattered-light spectra assuming e = 10~* (thin green-dashed), and e = 2 X 10~* (thin blue-dash-dotted), respectively. 
(C) and (D) show the normalized distributions of results from 4,059 bootstrap- re-samplings, for the data (thick red-solid), the simulated 
scattered-light spectra assuming e = 10"* {thin green- dashed), and e = 2 X 10~* (thin blue- dash- dotted), respectively. The three vertical 
lines mark the averages of the normalized distributions. AAf (i?) = is also shown as dotted lines for comparison. 



of the EW ratios between the scattered and stellar ab- 
sorption lines. We have made a comparison sample by 
Doppler-shift correcting the exposures with multiple sets 
of scrambled, instead of the real, redshifts. This com- 
parison sample has been passed through the same data 
pipeline, and provided the control comparisons, as well 
as the significance levels of any detections. Jackknife- 
and bootstrap- re-samplings have yielded consistent er- 
rors with those determined from the variations of the 
reference spectra. We have also produced spectra with 
simulated scattered signals, using the stellar absorption 
lines scaled by a factor of an assumed planet-to-star flux 
ratio. The simulated signals have been fully recovered 
through the data pipeline. Physical parameters of the 
system are required to have been determined, in order 
to predict the planet's relative radial velocity at each 
orbital phase, and further localize its imprinted signals 



on the recorded spectra. These parameters can be deter- 
mined from radial- velocity and photometric observations 
of transiting CEGPs with high precision and reliability. 
As a case study, we have analyzed very high S/N op- 
tical spectra of HD 209458b, acquired with the Subaru 
HDS in UT 2002 October. Our results suggest an av- 
erage planet-to-star flux ratio of (1.4±2.9)xl0~'*, in the 
wavelength from 554 nm to 681 nm, when the planet 
was during the orbital phases 11.0° < a < 33.9°. As- 
suming a Lambert-sphere phas e function (jPollack et al.l 
119861 : ICharbonneau et al.lll999f ). the equivalent geomet- 
ric albedo can be estimated as 0.8 ± 1.6, according to 
Equation [3l where we take the phase of the middle ex- 
posure (f> w 0.561, corresponding to a = 22.3°, out 
of 11.0° < a < 33.9° for all the exposures we used. 
The best previous existing constraint on scattered light 
from HD 209458b is based on direct imaging with the 
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MOST satellite ()Rowe et al.ll2006[ ). iRowe et all (|2006[ ) 
found a 1-cr limit on the planet-to-star flux ratio of 
4.88 X 10"*^, corresponding to a geometric albedo upper 
limit in MOST bandpass (400-700 nm) of 0.25. Our re- 
sults provide a useful point of comparison for this difficult 
measurement, covering a somewhat different wavelength 
range and based on different observational strategies and 
analysis techniques. 

These results do not seem yet able to support a com- 
parison to models of the HD 209458b atmosphere or to 
justify detailed physical interpretation, but they do not 
miss the mark by far and are thus quite encouraging. 
Specifically, the Subaru HDS data analyzed here consti- 
tute well under 5 hours of total exposure time; it would 
be possible to obtain more data, and there is no indica- 
tion in the present data of systematic problems (such as 
instrumental instabilities) that would prevent further in- 
tegration from producing a greater sensitivity. Moreover, 



one significantly more favorable transiting CEGP system 
(HD 189733) is known, and others may be discovered. In 
combination with measurements of re- emitted thermal 
radiation fe.g. lCharbonneau et al.|[2005D . scattered light 
studies offer the possibility of major qualitative advances 
in our understanding of CEGPs, a central riddle in the 
study of exoplanets generally. 



This work is based on data from the Subaru Telescope, 
which is operated by the National Astronomical Obser- 
vatory of Japan. We wish to recognize and acknowledge 
the very significant cultural role and reverence that the 
summit of Mauna Kea has always had within the indige- 
nous Hawaiian community. We are most fortunate to 
have the opportunity to conduct observations from this 
mountain. 



REFERENCES 



Brown, T. M., Charbonneau, D., Gilliland, R. L., Noyes, R. W., & 

Burrows, A. 2001, ApJ, 552, 699 
Burrows, A., Hubeny, I., & Sudarsky, D. 2005, ApJ, 625, L135 
Cameron, A. C, Horne, K., Penny, A., & James, D. 1999, Nature, 

402, 751 
Charbonneau, D., Allen, L. E., Megeath, S. T., Torres, C, Alonso, 

R., Brown, T. M., Gilliland, R. L., Latham, D. W., Mandushov, 

C, O'Donovan, F. T., & Sozzetti, A. 2005, ApJ, 626, 523 
Charbonneau, D., Brown, T. M., Latham, D. W., & Mayor, M. 

2000, ApJ, 529, L45 
Charbonneau, D., Noyes, R. W., Korzennik, S. C, Nisenson, P., 

Jha, S., Vogt, S. S., & Kibrick, R. L 1999, ApJ, 522, L145 
Chauvin, G. 2007, in Proceedings of the conference In the Spirit of 

Bernard Lyot: The Direct Detection of Planets and Circumstellar 

Disks in the 21st Century. June 04 - 08, 2007. University 

of California, Berkeley, CA, USA. Edited by Paul Kalas., ed. 

P. Kalas 
Collier Cameron, A., Horne, K., Penny, A., & Leigh, C. 2002, 

MNRAS, 330, 187 
Henry, G. W., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2000, 

ApJ, 529, L41 
Knutson, H. A., Charbonneau, D., Noyes, R. W., Brown, T. M., & 

Gilliland, R. L. 2007, ApJ, 655, 564 
Leigh, C, Cameron, A. C, & Guillot, T. 2003, MNRAS, 346, 890 
Marley, M. S., Gelino, C, Stephens, D., Lunine, J. I., & Freedman, 

R. 1999, ApJ, 513, 879 



Mayor, M. & Queloz, D. 1995, Nature, 378, 355 

Mazeh, T., Naef, D., Torres, G., Latham, D. W., Mayor, M., Beuzit, 
J.-L., Brown, T. M., Buchhave, L., Burnet, M., Carney, B. W., 
Charbonneau, D., Drukier, G. A., Laird, J. B., Pepe, F., Perrier, 
C, Queloz, D., Santos, N. C, Sivan, J.-P., Udry, S., & Zucker, 
S. 2000, ApJ, 532, L55 

Narita, N., Suto, Y., Winn, J. N., Turner, E. L., Aoki, W., Leigh, 

C. J., Sato, B., Tamura, M., & Yamada, T. 2005, PASJ, 57, 471 
Noguchi, K., Aoki, W., Kawanomoto, S., Ando, H., Honda, S., 

Izumiura, H., Kambe, E., Okita, K., Sadakane, K., Sato, B., 

Tajitsu, A., Takada-Hidai, T., Tanaka, W., Watanabe, E., & 

Yoshida, M. 2002, PASJ, 54, 855 
Ferryman, M. A. C. 2000, Reports of Progress in Physics, 63, 1209 
Pollack, J. B., Rages, K., Baines, K. H., Bergstralh, J. T., Wenkert, 

D., & Danielson, G. E. 1986, Icarus, 65, 442 
Rowe, J. F., Matthews, J. M., Seager, S., Kuschnig, R., Guenther, 

D. B., Moffat, A. F. J., Rucinski, S. M., Sasselov, D., Walker, 
G. A. H., & Weiss, W. W. 2006, ApJ, 646, 1241 

Seager, S. & Sasselov, D. D. 1998, ApJ, 502, L157-I- 
Seager, S., Whitney, B. A., & Sasselov, D. D. 2000, ApJ, 540, 504 
Sudarsky, D., Burrows, A., & Hubeny, I. 2003, ApJ, 588, 1121 
Sudarsky, D., Burrows, A., & Pinto, P. 2000, ApJ, 538, 885 
Winn, J. N., Suto, Y., Turner, E. L., Narita, N., Frye, B. L., Aoki, 
W., Sato, B., & Yamada, T. 2004, PASJ, 56, 655 



